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Введение. Исследуется процесс пиролиза пропана в про- 
точном реакторе в режиме подачи энергии через постоян- 
ный внешний нагрев. Для комплексного исследования 
процесса во всевозможных условиях и для обеспечения 
перехода от лабораторных установок к промышленным 
необходимо провести численный анализ результатов мо- 
делирования течения газа в реакторе с учётом химических 
процессов. В работе представлены результаты численного 
моделирования трехмерной динамики газового потока 
пиролиза пропана в реакторе в программном пакете 
АМЗУЗ Ешчеп с использованием компактной кинетиче- 
ской модели. 

Материалы и методы. Авторами предложена новая ком- 
пактная кинетическая схема пиролиза пропана, которая 
была получена на основе методов локального и глобаль- 
ного анализа чувствительности модели. Представлена 
математическая модель, используемая в пакете АМЗУЗ 
Ешепь которая представляет собой уравнения неразрыв- 
ности, сохранения импульса, сохранения энергии, допол- 
ненные уравнениями неразрывности для каждой компо- 
ненты газа. 

Результаты исследования. Впервые проведено численное 
моделирование трехмерной динамики газового потока пиро- 
лиза пропана в реакторе в программном пакете АМЗУЗ 
Ешеп с использованием компактной кинетической модели. 
Проведены расчеты динамики газового потока пиролиза 
пропана в лабораторном реакторе с учетом процессов диф- 
фузии, химических реакций и их тепловых эффектов. 
Результаты численных расчетов и экспериментальные ис- 
следования по конверсии пропана хорошо согласуются 
между собой. 

Обсуждение и заключения. Результаты исследования и 
моделирования пиролиза пропана могут составить основу 
для описания процесса в объеме реактора под воздей- 
ствием лазерного излучения. 
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Введение. В Институте катализа им. Г.К. Борескова СО РАН с целью создания экономически эффек- 
тивной технологии переработки углеводородного сырья изучают химическую реакцию пиролиза пропана в ла- 
бораторном реакторе проточного типа. Газофазные условия при протекании процесса формируются за счет 
внешнего нагрева стенок реактора и за счет подачи лазерного излучения, что обеспечивает дополнительную 
генерацию радикалов [1]. При этом происходит повышение реакционной способности смеси, и тем самым уве- 
личивается выход полезных продуктов. Кроме того, уменьшается образование кокса и побочных продуктов 
реакции вследствие возможности проведения процесса пиролиза пропана при более низкой температуре. 

Для комплексного исследования процесса при различных условиях и для обеспечения перехода от ла- 
бораторных установок к промышленным, необходимо провести численный анализ результатов моделирования 
течения газа в реакторе с учётом химических процессов. Общемировой практикой при математическом моде- 
лировании является использование мощных программных комплексов ЗО-моделирования промышленных и 
лабораторных реакторов (например, коммерческие пакеты Еепё, Е]о\1тага, Е№юо\/У151оп), которые могут экс- 
плуатироваться на однопроцессорных и многопроцессорных ЭВМ. 

При трехмерном моделировании газовой динамики в реакторе основную трудность представляет собой 
выбор кинетической модели реакции. Пиролиз пропана является радикально-цепной реакцией, и численное 
интегрирование получаемой жесткой системы дифференциальных уравнений представляет собой вычислитель- 
но затратную задачу. Для описания пиролиза пропана были установлены многочисленные кинетические моде- 
ли, среди которых можно условно выделить компактные (описывают пиролиз пропана в заданных условиях с 
относительно небольшим числом стадий, до 20) и детализированные. Для ЗО-моделирования, как правило, ис- 
пользуют компактные схемы, что обусловлено большими временными затратами, возникающими при расчетах 
по детализированным схемам [2, 3]. 

В настоящее время в основном изучено одно- и двумерное моделирование процесса газофазного пиро- 
лиза пропана [4—6]. Для численного изучения процессов теплообмена и динамики данного процесса в трехмер- 
ном пространстве авторами был выбран пакет АМЗУЗ Ещепй (лицензия Сибирского суперкомпьютерного цен- 
тра СО РАН) [7]. 

В настоящей статье рассматриваются результаты верификации численной модели трехмерного реаги- 
рующего газа, построенной на основе пакета АМЗУЗ Ешеп путем сравнения с экспериментальными данными, 
полученными для термического пиролиза пропана в проточном реакторе в режиме подачи энергии в реактор 
через постоянный внешний нагрев. 

Кинетическая модель. Термическое разложение предельных углеводородов сопровождается разнооб- 
разными и многочисленными параллельными процессами, что приводит к образованию десятка веществ и воз- 
можности протекания сотни реакций. Количество стадий и веществ, включенных в кинетическую модель, 
должно быть ограничено. Это мотивируется необходимостью интегрирования кинетической модели в про- 
граммный комплекс ЗО-моделирования АМЗУЗ Ешеп{ для расчета трехмерной динамики газа. Брутто-схемы не 
дают удовлетворительного описания пространственного распределения по реактору протекания химических 
реакций с радикально-цепными механизмами. Использование детальных кинетических механизмов при моде- 
лировании процесса пиролиза пропана, включающих химические превращения, важно для полного и точного 
описания процессов в широком диапазоне температур. Однако их использование в трехмерных численных рас- 
четах реакторов приводит к неприемлемым временам расчетов даже на современных суперкомпьютерах. Таким 
образом необходимо найти баланс между точностью модели и временем вычислений. 

Рассмотрев различные кинетические схемы в пакете АМЗУ$З, была получена схема приемлемого раз- 
мера(табл. 1), описывающая расход пропана и образования основных продуктов реакции [8—11]. Константы 
элементарных стадий широко представлены в справочной литературе, в научных статьях, а также в химических 
базах данных, прежде всего в МТЗТ. В таблице 1 приведены элементарные химические стадии пиролиза пропа- 
на. 
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Таблица 1 
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Математическая модель. Математическая модель представляет собой уравнения неразрывности, со- 


хранения импульса, сохранения энергии, дополненные уравнениями неразрывности для каждой компоненты 
газа: 
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Уравнение состояния смеси идеальных газов имеет вид: 
р-РАТХ О, М.) 


Начальные условия задают состояние рассматриваемого объекта в начальный момент времени: 


р(х,&, ) = рь, У(х,к) = (и(х&),у(хк),м(хь))", Рок) = р, Г(хь)=Т, Г(хь)=У. 





3 > 
Здесь р — плотность газа (кг/м”), У — скорость газа (м/с), р — давление (Па), рЁЕ’ — полная энергия в 
единице объема (Дж/м3), Й; — энтальпия (Дж/кг), с, — удельная теплоемкость вещества / (Дж/Кэкг), 
п, 3 
Е К, — нагрев (охлаждение) в результате химической реакции и внешних воздействий (Дж/сэм”), 
м 
и 
С № ть . 
У, — локальная массовая доля вещества 1, К, =М„>К, — объемная скорость образования вещества 1 
и=1 


(кг/сэм”), Е, — молярный коэффициент Аррениуса образования (расхода) вещества { в ходе реакции 7 (1/сэм°), 


М ил — Молекулярная масса вещества 1 (кг), ^ — универсальная газовая постоянная (Дж/мольэК), Г — темпе- 


ратура (К). 
Более подробно математическая модель представлена в работе [7]. 


Результаты исследования. Были проведены расчеты динамики газового потока пиролиза пропана с 
учетом процессов диффузии, тепловых эффектов реакции и тепловых процессов с включением разработанной 
кинетической модели. Кинетическая модель является первостепенной задачей при моделировании любого 
нефтехимического процесса. На ее основе получают оптимальные условия проведения нефтехимического про- 
цесса. Исследуются важные явления, например, реакционная способность участвующих в реакции всех ве- 
ществ, образование и устойчивое существование индукционных (взрывных) периодов и т. д. [12, 13]. В матема- 
тической модели реактора рассматривается нестационарное пространственное трехмерное течение газа; лами- 
нарный режим течения с диффузией процессов. В модель включены уравнения сохранения массы, импульса, 
энергии и уравнения сохранения отдельных химических веществ. Скорость реакции рассчитывается из уравне- 
ния Аррениуса на основе кинетической модели, которая был разработана при использовании методов локаль- 
ного и глобального анализа чувствительности [ 14—17]. 

Температура смеси, подаваемой в реактор, составляла 600 К. Пристеночная температура задавалась в 
диапазоне от 820 до 990 К, температура входных и выходных трубок — 300 К. Подаваемые расходы через 
верхний и нижний входы — 1,758 мг/сек. Ввод защитного газа (метана) осуществлялся с расходом 0,221 мг/сек 
и температурой 300 К. 

Результаты моделирования (профиль температуры и мольных долей основных компонентов) приведе- 
ны на рисунках 1—5 для одного из моментов установления течения в реакторе, когда распределение температу- 
ры и скорости потока реагентов на выходе реактора меняются незначительно. 

Защитный газ (метан) ограничивает область реакционной зоны и предохраняет от перегрева торцевые 
стенки с окнами для излучения. Метан преимущественно остается в буферных зонах, и лишь малая его часть 
проникает в реакционную зону. Потоки холодного газа, поступающего через вводы исходной смеси, не успе- 
вают охладить двигающуюся по реакционной зоне смесь. Температура газа в реакционной зоне меняется от 
высоких пристеночных значений до минимума в центре. Такое распределение способствует потоку энергии от 
стенок к оси реактора. Химическая реакция протекает почти по всему объему реакционной зоны с более актив- 
ным превращением в пристеночной разогретой области. Появление пропана и продуктов реакции в правой и 
левой «защитных» областях реактора объясняется диффузией газа. Большая массовая доля продуктов реакции в 
правой части реакционной зоны объясняется тем, что в продолжающем движение газе реакции не прекращают- 
ся. Превращение пропана проходит во всей реакционной зоне. 


И. М. Губайдуллин и др. Численный анализ динамики газового потока пиролиза пропана 
1. М. СиБау4иШт апа йе ошйегх. Митетса[! апа[р?515 ор газ По 4упатгсу орргорапе руго[у515 








3002402 4.01е+02 5.02е+02 503е+02 7.042+02 8.05е+02 9.0Бе+02 9,73е+02 


Рис. 1. Распределение температуры (К) в реакторе в продольном сечении 
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Рис. 2. Распределение в продольном сечении мольных долей пропана 
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Рис. 3. Распределение в продольном сечении мольных долей метана 
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Рис. 4. Распределение в продольном сечении мольных долей этилена 
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Рис. 5. Распределение в продольном сечении мольных долей водорода 
Ето. 5. Нуаговеп то[е гасйоп @триноп а опепи4 та! 5еспоп 


Для верификации трехмерной модели использовалось сопоставление экспериментальных и расчетных 
данных по зависимости конверсии пропана от пристеночной температуры. Конверсия пропана при таком рас- 
пределении температуры и с приведенной выше схемой близка к экспериментальным значениям (табл. 2). 


Таблица 2 
ТаБЕ 2 


Сравнение расчетных и экспериментальных данных по конверсии пропана 


Сотрапбоп оГсасшеиеа апа ехреттепк ааа оп ргорапе сопуегяюоп 
Т К 
Тетрегайтге, К 
0 
р о 73 127 22 36 
Ехрегитепь % 
Расчет, % 
И 6,8 15,6 В. 37,1 
Са|сШайноп, % 


Заключение. В программном пакете АМЗУЗ Еег{ проведены расчеты динамики газового потока пи- 
ролиза пропана в лабораторном реакторе с учетом процессов диффузии и тепловых эффектов реакции в про- 
граммном пакете АМЗУ5З Ееп{ с включением разработанной компактной схемы реакции. Проведен вычисли- 
тельный эксперимент и найдены зависимости выхода продуктов реакции и расхода подаваемой смеси по длине 
реакционной зоны. Результаты численных расчетов и экспериментальные исследования по конверсии пропана 
хорошо согласуются между собой. Успешные результаты моделирования неизотермического химического ре- 
актора с ламинарным потоком реагентов открывают возможность использования этого типа реакторов для изу- 





чения процесса в объеме реактора под воздействием лазерного излучения. 
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